oo 

O ■ Evaluation of the semiclassical coherent state 
(N . . 

qH propagator in the presence of phase space caustics 

A D Ribeiro^'^ and M A M de Aguiar^ 

^ Institute de Fisica, Universidade de Sao Paulo, Usp, 05315-970, Sao Paulo, Sao Paulo, Brazil 
^ Institute de Fisica "Gleb Wataghin", Universidade Estadual de Campinas, Unicamp, 
, 13083-970, Campinas, Sao Paulo, Brazil 

E-mail: aribeiroOif i .unicamp .br ,aguiar(3ifi .unicamp.br 

> 

. Abstract. A uniform approximation for the coherent state propagator, valid in the vicinity 

^ ' of phase space caustics, was recently obtained using the Maslov method combined with a dual 

representation for coherent states. In this paper we review the derivation of this formula and 
apply it to two model systems: the one-dimensional quartic oscillator and a two-dimensional 
chaotic system. 
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1. Introduction 

The representation of coherent states has been used to describe a wide variety of physical 
CD I systems [1]. In particular, coherent states provide a natural phase space representation of 

quantum mechanics and are specially well suited to the study of the semiclassical limit. The 
coherent state representation was first formalized by Bargmann in 1961 [2j and later used by 
Glauber [3] to describe the electromagnetic field in quantum electrodynamics. Fairly complete 
^ ■ review articles on coherent states and applications can be found in references [Hill [5]. 

■ The quantum propagator K{z"* , z' , T) = {z"\e~^^'^/^\z') represents the probability amplitude 

that an initial coherent state \z') evolves into another coherent state \z") after a time T. A path 
integral formulation for this propagator was introduced by Klauder [6] . The paths contributing to 
K{z"*,z',T) are those connecting {q' ,p') = {{z'\q\z'), {z'\p\z')) to {q",p") = {{z"\q\z"), {z"\p\z")). 
In the semiclassical limit, it turns out that the most important paths are complex classical 
trajectories governed by the hamiltonian function H = {z\H\z), with boundary conditions 
involving the average values q' , p' , q" and p" . Klauder was the first to consider this type of 
approximation [7], being followed by a number of other contributors [51 El |TD]. More recently, 
a detailed derivation of the semiclassical coherent state propagator for one dimensional systems 
was published [TT] . In the last two decades much numerical and analytical work has been done in 
semiclassical methods with coherent states for one [l2l [111 H [ig m HZl [H US] and two pnll2T] 
dimensional systems. A recent review can be found in reference |22j . 

While the semiclassical propagator is usually very accurate for short times, phase space 
caustics and non-contributing trajectories inevitably appear as T increases, introducing large 
errors and imperfections [121 US (Ml EHl El]. Non-contributing trajectories must be identified 
and excluded from the calculation because their contributions to the propagator are non- 
physical. From the mathematical point of view, non-contributing trajectories correspond to 
forbidden deformations of the contour of integration necessary to carry out the stationary 
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phase approximation that leads to the semiclassical formula. Phase space caustics, on the 
other hand, are special points where the amplitude of the semiclassical propagator diverges and 
the approximation simply breaks down. In order to calculate the propagator in the vicinity of 
caustics one needs to improve the semiclassical approximation, going beyond the usual quadratic 
expansion. 

In two recent papers [231 124j we have derived a uniform approximation for the coherent state 
propagator that remains finite in the presence of phase space caustics. The derivation involved 
the introduction of a dual representation for the coherent states and the method of Maslov [25j. 
In the present article, we review the formalism used in these previous papers and apply it to 
the one-dimensional quartic oscillator and to the two dimensional chaotic Nelson potential. We 
show that the uniform formula completely eliminates the divergences caused by the caustics, 
providing a very accurate semiclassical description of the propagator in these regions. 

This paper is organized as follows: in the next section we briefly review the representation of 
coherent states and the quantum propagator. Section[3]describes the semiclassical approximation 
to the propagator based on a second order expansion around stationary trajectories. The dual 
representation for coherent states is introduced in section U] and used in section [5] to derive the 
uniform approximation. In sections [6] and [7] we present numerical applications of the uniform 
formula and, in section [8l we present our final remarks. 



2. The coherent state propagator 

Let Hq = Tiijj{ci^a + 1/2) be the Hamiltonian of a harmonic oscillator of mass m and frequency 
oj. The normalized coherent states of Hq are defined by [H [U [5] 

\z) =e-^\'\\'^^\Q), (1) 



where 

1 f q .p\ ^ f 1 -P 



and |0) is the oscillator's ground state. The real labels q and p are the average values 
of the position and momentum operators respectively and the length scales b = y/h/[muj) 
and c = V mhuj satisfy be = h. Three important properties of the coherent states are 
(over)completeness, overlap relation and eigenvalue equation: 

I ^d'^z f .dqdp 

k)— (^i=y (3) 

{zi\zj) = expj-^lzip + z*Zj - ^kjf } (4) 

and 

a\z) = z\z). (5) 

It will also be important in the derivation of our uniform approximation to define non- 
normalized coherent states, or Bargmann states [2j, by 

k) = e^''|0>. (6) 



For these states the unit operator and the overlap equation become 



\z) {z\dz and = e^* (7) 



The quantum propagator in the Bargmann and the coherent states representations are given, 
respectively, by 

k{z"*,z',T) = {z"\e-'^^/^\z') (8) 

and 

K{z"*,z',T) = {z"\e-'''^/V) = e--2\''\"~^^\'"\\{z"\z',T). (9) 

These two quantities contain the same physical information and differ only in the normalization. 
For the purposes of the theory to be developed in sections H] and [5] it shall be more convenient 
to work with the Bargmann states. 

The propagator K(z"* , z' , T) can be written in terms of path integrals, from which standard 
semiclassical approximations can be performed. Here we present a very brief summary of path 
integral formulation, referring to [llj for the details. The first step is to divide the propagation 
time T into N small intervals of size e = T/N so that 

K{z"\z\T)= lim {z"\ e-'^'l^ . . .e-'^'/^, \z'). (10) 

N^oo V ' 

A'' times 

Next, the coherent state unity operator ([4]) is inserted between each infinitesimal operator 
^-iHe/h^ The path integral formula, obtained by evaluating the expression for each resulting 
infinitesimal propagator {zjj^i\e~'^^'^/^\zj) , reads as 



K{z"*,z',T)= lim / TT l^ L«^'=-» in ^ ^ 



N^oo 




where we have defined i/fc_j_i/2 = (•ZA;+i|-f^|-Zfc)/(-2fc+i|-Zfc) and identified \z') = \zo) and {z"\ = {zj\f\. 
Equation (llip represents the path integral formula of the quantum propagator. Written as a 
function of the numbers {qj,Pj), the propagator becomes an infinite sum of contributions of all 
possible phase space paths linking the initial point {q' ,p') to the final {q",p")- 

3. Second order semiclassical approximation 

By taking the formal semiclassical limit fi ^ 0, the integrals (|lip can be performed [TT] by the 
saddle point method [26j. It can be shown that critical paths, those whose contribution to the 
integral are relevant, are classical trajectories {Q{t), P{t)), satisfying the boundary conditions 

b c b c b c b c 

and governed by the average hamiltonian H = {z\H\z). One might think initially that the 
trajectory starting at {Q{0), P{0)) = {q',p') and ending at {Q{T), P{T)) = {q",p") would be the 
only solution to these equations. However, these boundary conditions are very restrictive and 
such a trajectory usually does not exist: indeed, giving the initial position and initial momentum, 
the trajectory is completely determined so that the final point {Q(T), P{T)) is generally different 
from {q",p"). This means that, in general, there is no real critical path to the integral (jlip . 
Complex trajectories, however, can usually be found if we analytically extend the integration to 
the complex phase space, letting Q and P be complex variables. 

In this case it is more convenient to introduce new variables u and v, instead of using the 
complex position Q and momentum P, such that 

u=—;=\—+i—] and V = —;= [ «— • (13) 



In terms of u and v the classical equations of motion become 



u 



1 dH 

ih dv 



— —— and V 



1 dH 

ih du ' 



(14) 



where H = {v\H\u), and the boundary conditions assume a simpler form, 

uiO): 

v{T) 



j_ /Q(r) _ .p(T)\ ^J_ (<f_ 

y/2\ h ' C ) 



.P" 



(15) 



Notice that this does not imply that v{Q) 



and u{T) = z" . Given n(0) = z' and v{T) 



a complex trajectory {u{t),v{t)) can be calculated and the values of u(0) and u{T) come out 
of this calculation. In general there might be more than one trajectory governed by ()14p and 
satisfying ([15]). 

Returning to Eq. ([TT]) . by expanding the exponent up to second order around the complex 
classical path, we find the following semiclassical formula for the propagator 



KfXz"\z\T) 




exp 



{i(5 + e)-i(|.? + |.T)} 



(16) 



where the index (2) means "second order expansion" . The sum in Eq. ()16|) is over the complex 
classical trajectories, as discussed, and 
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Finally M^y is an element of the tangent matrix M defined by 



6u{T) 
5v{T) 



Muu Muv 



<5n(0) 



(17) 
(18) 

(19) 



where 5u and 5v are small displacements around the classical trajectory. The elements of the 
tangent matrix are related to second derivatives of the action S [11]. The phase of M„„ contains 
important Maslov phases. 

The corresponding semiclassical formula for k{z"* , z' ,T) is identical, except that it does not 
have the normalization factor in the exponent. 

In the strict limit where h ^ 0, the semiclassical propagator becomes a delta function at 
the phase space point {q",p") linked to {q',p') by a real classical trajectory. Therefore, for 
small but finite h, we expect large contributions to the propagator arising from nearly real 
trajectories. The more the trajectory wanders into the complex p and q space, the less it 
should contribute to Ksc (z"* , z' ,T). By inspection of Eq. (fT6j) . we see that this statement is 
true only if the real part of the exponent in ()16p is negative. Trajectories for which such real 
part is positive furnish non-physical contributions to the propagator that become arbitrarily 
large when h goes to zero. These are non-contributing trajectories and must be excluded from 
the calculation. They correspond to complex critical paths whose steepest descent contour of 
integration cannot be reached by deformations allowed by Cauchy's theorem. Their structure 
are closely related to the Stokes Phenomenon [27\ I28j . discussed in a number of papers on 



semiclassical approximations [121 [131 [201 [2T1 [291 [301 [3T] . Iii particular, Ref. [31] shows an explicit 
example where non-contributing trajectories arise from forbidden deformations of the original 
contour of integration. 

Another common source of complications in semiclassical formulas are focal points or 
caustics. These are special points where the second order approximation breaks down because 
of singularities in the formula's pre-factor. In the case of Ksc {z"* , z' ,T) this happens when 
-/Vft,„ = 0. According to Eq. (jl9p . if goes to 0, we can set small initial displacements 
5u(fi) = and 5v(Q) ^ such that 6u{T) / and 6v{T) = 0, implying that there are at least 
two nearby trajectories satisfying the correct boundary conditions (jlSp . The point where these 
trajectories coalesce is called a phase space caustic and the semiclassical formula (jl6p fails in its 
vicinity. Contrary to non-contributing trajectories, trajectories going through caustics cannot 
simply be excluded, since the problem lies on the approach used, and not on the orbit itself. 
Thus, to evaluate the semiclassical propagator in the vicinity of phase space caustics one needs 
better approximations, beyond the second order. We derive such an approximation in the next 
two sections. 



4. Dual representation for the coherent state propagator 

The most direct way to improve the quadratic approximation is to go back to Eq. (jllh and expand 
the exponent of the propagator to third order around the critical paths. This, however, would 
be extremely complicated, since the integral (|lip is multi-dimensional. A simpler solution is to 
follow the method proposed by Maslov [25]. To illustrate the idea, suppose that a semiclassical 
approximation for ^lJ{q) has a singularity at q = qq. If we know the corresponding semiclassical 
formula for in the momentum representation we can write 

i^sM = J {qo\p)Mp)dp = J ^sc{jp)e'"^''^dp. (20) 

If the integral over p is performed by the usual second order stationary phase approximation, the 
singularity in tpsdQo) is recovered. However, doing a third order stationary phase approximation 
produces a more accurate expression involving an Airy function which remains finite at qq [32j . 
Therefore, Maslov's method consists in finding the desired semiclassical approximation in a 
conjugate representation and transforming back to the original one by a third order expansion. 

The problem in applying this idea to coherent states is that they are defined in phase space 
and do not have a natural dual representation. Since z and z* play the role of conjugate variables, 
we can think of z as g and z* as p, and the coherent state propagator is always written in the 
mixed p — q representation. It is impossible to write it in q — q or p — p forms, since one cannot 
write a matrix element with two kets or two bras. To overcome this difficulty we need to define 
a proper linear application to play the role of the dual representation and we shall do that using 
the non-normalized Bargmann representation. 

Given the propagator k{z"* , z' ,T), the associated dual propagator is defined as 

~k{z,z',T) = ^= [ k{z"*,z',T)e-'""'dz"*, (21) 
V27rz Jc 

and the inverse application by 

k{z"*,z',T) = ^= I k{z,z',T)e'"*'dz, 
V27ri Jc 

where C and C are convenient paths as specified in Ref. [23|. These equations are reminiscent 
of (|20p . Notice that equation (j2ip can also be written as 

k{z, z ,T) = —= I — — dz , 

\/2Tn Jc {z \z) 



(22) 



which can be interpreted as an attempt to 'cancel' the bra {z"\ and introduce another ket \z). Of 
course k is not a matrix element and, therefore, the above application is not a true representation. 

Equation ()22p is the starting point to make improvements on ksc {z"* , z' ,T). In the regions 
where both propagators are free of caustics, the semiclassical version of A;(2;, z' ,T) can be obtained 
by performing the integral (j2ip using the standard second order saddle point method with 

k{z"*,z',T) replaced by kil\z"* , z' ,T). The result is [23] 



{z, z', T) = Y.J^ exp l'-S{z, z',T) + ^g{z, z' ,T)\ . (23) 

traj. V ™ 

The trajectories summed in this equation are not the same as those of Eq. (|16p . since they satisfy 
the conditions u(0) = z' and u{T) = z. As usual, ^(0) and v{T) are not fixed, but come out 
of the integration of Hamilton's equations with the above boundary conditions. The tangent 
matrix element M„„ is given by Eq. (1190 and g{z,z',T) is the function Q calculated with the 
new trajectory. The new action S{z, z' ,T) = S{z"* , z' , T) + ifizz"* is the Legendre transform of 
5, where z"* = z"* {z, z' ,T) is obtained from the relation —ihz = dS/dz"* . 

According to Eq. (jl9p . when M^^ is zero, M^v is usually not zero, which implies that if 
ki^ {z"* , z' ,T) has a caustic for a complex trajectory satisfying n(0) = z', v{T) = z"* , implying 

values of u{T) and i'(O), then ksc {z, z' ,T) will not have a caustic when calculated at the same 
trajectory, i.e., for u{0) = z' and z = u{T). 
Inserting ([23]) into ([22]) leads to 



ksciz"*,z', T) = ^i J-^ ,iSi.,.'^H^g+.."^.,.'^) (24) 
V27rz Jc V ^uv 

where we omit the sum for simplicity. Eq. (j24p is an integral representation for the semiclassical 
coherent state propagator. In order to calculate it we need to sum the contributions of all 
trajectories starting at m(0) = z' and ending at u{T) = z lying on the curve C. Eq. (fTUj) is 
recovered if the second order saddle point method is applied to (fM|) . the critical paths being 
precisely the orbits given by Eqs. (fT^ and (fTSj) . 

Expanding the exponent of Eq. (I24p to third order around the critical paths leads to the 
so called regular approximation derived in [24J. The regular formula provides satisfactory 
results only if the critical trajectories are not too close to caustics, so that each trajectory still 
contributes independently to the propagator. For the transitional approximation, the exponent 
of (|24p is expanded around the trajectory that lies exactly at the phase space caustic. Since this 
trajectory is not a critical one, the formula is good only if critical trajectories are sufficiently 
close to the caustic. In this paper we will be concerned only with uniform approximations, which 
provide a global semiclassical formula, reasonably accurate over all the space spanned by the 
parameters z' , z"* and T. 

5. Uniform Approximation 

A uniform approximation for the coherent state propagator was obtained in Ref. [23j for one- 
dimensional systems and in Ref. [21] for two dimensions. The expression presented here is 
slightly different from those in Refs. |23l I24j . as we point out below. 

The simplest type of singularity that can appear in the semiclassical propagator occurs 
when two nearby trajectories coalesce at the caustic. In this case the function (p{z) = 
i{S{z, z' ,T) — ihzz"*)/h has two stationary points (corresponding to two complex trajectories 
satisfying the same boundary conditions) that coalesce as z"* (here considered as a parameter) 
approaches the caustic. The basic idea of the uniform approximation is to map this complicated 



function into a simpler function with the same critical points and the same behavior in the 
neighborhoods of these points [33]. Since all that matters in the semiclassical Hmit is the 
neighborhood of the critical points, the integral with the new function should give about the 
same results as the integral with the original function. 

In the case of two coalescing trajectories the appropriate function is N{x) = A — Bx + ^x^, 
which has two saddle points at ib\/B that coalesce as ;B ^ 0. Therefore we write [see Eq. (j24p ] 

V{z,z',T)ei^'-'''''^^+''"' dz = J J{x)e^^^Ux, (25) 

where the function J{x) includes the Jacobian of the transformation z ^ x and the contribution 

of the smooth term V{z,z',T) = ''^e^^^^'^''^^. In Refs. [Ml [Ml the logarithm of this term 
was also included in the function (p{z). However, since V varies slowly with h it is reasonable to 
leave it out. From the numerical point of view it turns out that the present prescription is also 
more accurate than the ones in |23i . 

Imposing that the value of N{x) at the saddle points ib\/F coincide with the value of (f> at 
critical points 2:1^2 (related to the two critical trajectories of dM])) we find 



A = —{Si+S2) and B 



2/3 



(26) 



where Si is the complex action of the trajectory related to Zi. This determines completely the 
function N(x). 

Next we impose the equivalence between N{x), in the vicinity of x = =b\/^, and 4>{z), in the 
vicinity of z = z\^2- This is done by writing 



dx"^ 



5x^ = ±2VB 5x^ 



which, by identifying 



dz"^ 



bz^, 



21,2 



dz 



ih 



21,2 



(27) 



(28) 



21,2 



implies that the jacobian of the transformation calculated at each saddle point, ji^2; amounts to 

" \ ih (M,,/M„,)|,^/ ^^^^ 
The function J{x), therefore, can be conveniently written as 



J{^) = ^ (-^1 + J2) - ^ (-^2 - Jl) + 0{x'), 



(30) 



where J,- 



Muv ' e 



/2 ±g(.-,.',T) 



Finally, we write the uniform formula for the propagator as 

k^^{z"\z',T) = f J{x) e-^-^"+"'/=^d: 

^/2■Ki J 

or, discarding global phases. 



ji, so that when x xi^2 = ±V^, we have J{xi^2) = Ji,2- 

(31) 



K^{z"\z\T) 



vr e 



A 
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Vb 



F[{B) + {g,+g2)F,{B)\, 



(32) 



where gi,2 = ^J^VB/ {M^. 
and Fi is given by 



, F( is the derivative of Fi with respect to its argument 



^i{W) = — f dtexp 



1 



Wt + -t^ 
3 



(33) 



for i = 1,2,3. The index i refers to three possible paths of integration Cj, giving rise to three 
different Airy functions [26j. The correct path is determined by Cauchy's theorem, but, in 
practice, we can use physical criteria to justify the choice of Q. The normalized propagator 



K^^{z"*,z',T) is obtained multiplying k^^ by e 



To finish this section we emphasize that the application of the uniform semiclassical formula 
(I32p always involves the two complex trajectories coalescing at the caustic and the choice of 
a contour of integration. In the next two sections we present numerical examples where these 
trajectories and the integration path will be pointed out explicitly. 



6. The quartic oscillator 

As a first application of the uniform approximation we consider the one-dimensional quartic 
oscillator. This simple system illustrates well the problems of non-contributing trajectories and 
phase space caustics. Moreover, the structure of the complex phase space is easy to visualize and 
shows clearly the complications that arise as the propagation time increases. The hamiltonian 
is 

^ = 1^'' + + (34) 

and we set B = 0.1 and % = 1.0. The initial state \z') is chosen at q' = 0, p' = —2.0 with 
b = 1.0, while the final \z") is given by q" = 0.5, p" = 0.5 with the same width b = 1.0. The 

calculation of {z"* , z' , T) for these fixed states as function of T shows that, for T ~ 2.5, the 
semiclassical result has an unphysical peak, revealing the presence of a caustic. We will show 
that this inaccurate result can be controlled with the uniform formula ()32p . 

In order to find complex trajectories satisfying the boundary conditions (115p we follow the 
recipe of Rubin and Klauder [13] and define 

g(0) =q' + w and P(0) = p' + i-w, (35) 

c 

where w = a + ifi is a complex number. It is easy to check that this choice of initial conditions 
satisfy the first of the boundary conditions ()15p for all values of w. The idea, therefore, is to 
propagate trajectories for all possible w, picking those satisfying the second of the boundary 
conditions (fT5]) . Specifically, for each w we propagate the complex trajectory starting from ([35]) 
and calculate 

1 (QiT) P{T)\_ 1 fQ" P"\ 

where Q{T) and P{T) are complex and Q" and P" are real variables, obtained by taking the real 
and imaginary parts of vt- The trajectories for which {Q",P") = {q",p") are the ones needed 

to calculate Kil\z"* , z' ,T) and K^^ {z"* , z' ,T) . 

Notice that the origin a = /3 = corresponds to the real trajectory starting from (5(0) = q' ■, 
P(0) = p' . Therefore, the larger the \vj\ the more complex is the corresponding trajectory and 
the smaller its contribution to the propagator. Thus, we can restrict our search to a small 
vicinity of the origin in the complex w plane. 

By fixing , p' and T, the values of the resulting Q" and P" can be seen as function of w. 
In Fig. [T] we represent the curves of constant Q" superimposed with those of constant P" for 
different values of T in the (a, /?) plane. Stars identify the points where Q" = ({' = 0.5 and 





Figure 1. Curves of constant Q" superimposed with those of constant P" in the w plane. 
Dashed hnes refer to Q" = and the dotted hnes to P" = 0. From panel (a) to (f), the value 
of T is to T = 0.06, 0.24, 0.70, 1.02, 2.20 and 2.70, respectively. Stars are centered at the point 
Q" = P" = 0.5. White stars represent trajectories belonging to the family /I, grey stars to 
/2, and the black stars to /3. The circles in panel (e) also represent potential candidates to 
be included in the sum (|16p . illustrating the existence of other trajectories satisfying the proper 
boundary conditions (jlSp . In this case, since they lie far from the origin, they were not included 
in the calculation. 
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Figure 2. Exact propagator (full line in all panels) and semiclassical results for Ksc {z"* , z' ,T) 
and Kg^{z"* , z' ,T) for the quartic oscillator. Panel (a) shows the exact result and the individual 

(2) 

contribution of each family of complex trajectories to Ksc ■ /I (dashed), /2 (dotted) and /3 
(dash-dotted). Panel (b) shows the propagator for T < 1.8 and displays the exact result, the 
contributions of /I (dashed) and /2 (dotted) and the combined contributions of /I and /2 
(dash-dotted). The uniform propagator iC"" calculated with /I and /2 is also shown, evaluated 
for each possible path: Ci (squares), C2 (circles) and C3 (triangles). Panel (c) focuses on the 
interval 1.8 < T < 4.2 and shows the exact result, the contributions of /2 (dotted) and /3 
(dash-dotted) and the combined contributions of /2 and /3 (dashed). The uniform propagator 
calculated with /2 and /3 is also shown along the three independent paths. Panel (d) shows 
the final result obtained with the uniform formula (dots) and exact result. 

P" = p" = 0.5, our desired final values. As T is changed, the position of the stars move in the 
w plane, forming families of contributing trajectories. Circles (only in Fig. [I{e)) are the same 
as stars, but their contributions were not included because they lie very far from the origin. 

White stars in Fig. [1] represent trajectories belonging to the family /I, which are close to the 
origin when T = but move away as T increases. Grey stars refer to the family /2. Trajectories 
belonging to this family start off quite complex, but approach the origin of the w plane as 
T tends to 3.2, approximately. After that, they move away from the origin again. Finally, 
trajectories of the family /3 are represented by black stars. They give important contributions 



for T > 2. Other trajectories satisfying these same boundary conditions exist (see Fig. HJe)) but 
their contribution to the propagator is not significant. 

For T = 0, a contour plot hke those in Fig. [1] would display a grid of straight lines, vertical 
for Q" = const, and horizontal for P" = const., with the only contributing trajectory lying at 
a = 1/4, P = —5/4. For short values of T, as in Fig. [T][|a) for T = 0.06, the grid deforms slightly 
and the contributing trajectory moves away from the origin. For T = 0.24, Fig. [IJb), there is 
still a single contributing trajectory, but two defects on the grid are seen approaching the origin. 
As discussed in Refs. [131 SO]) these defects are critical points of the map vt = vt{w), where 
~ ~ ^' addition, if the second derivatives of Q" and P" are non-zero, it can be shown 
that the map becomes two-to-one in the vicinity of the defects, implying two different trajectories 
(corresponding to two distinct initial conditions) satisfying the same boundary conditions (jlSp . 
These defects, therefore, can be identified with the phase space caustics. 

For T = 0.70 and T = 1.02, panels (c) and (d), two complex trajectories can be seen and for 
T = 2.20 and T = 2.70, panels (e) and (f), three of them are close to the origin. For T = 2.70, 
the trajectories belonging to /2 and /3 are the closest they get to each other, moving away for 
larger T. The effect of the nearby caustic should be pronounced at this point. 

Figure [2]J^a) shows the exact propagator and the individual contribution of each family, /I, /2 

(2) 

and /3, to the propag ator Krc{z"\z',T) for < T < 6. For very short times the contribution 
of /I is clearly dominant, and reproduces very well the exact results by itself. For T > 0.5 
the family /2 becomes more and more important, whereas the contribution of /I decreases. 
However, for T k, 1.0 the exact curve is not reproduced by /I nor by /2, suggesting that both 
families should be included simultaneously. 

FigureEJb) shows a detailed plot for the interval < T < 1.8 of such a combined contribution. 
The result obtained is clearly not accurate. However, by looking at Fig. [H we notice that these 
trajectories are close to a phase space caustic, and, therefore, a divergent behavior is expected. 
Results for Kg^{z"* , z' ,T) are presented in the same panel (b) for each possible integration 
path, Ci, C2 and C3 (see equation ([33|) ). Choosing the path Ci for < T < 0.65, and C3 for 
0.65 <T< 1.8, the exact result is satisfactorily reproduced. 

Returning to Fig. [2ja), we noticed that, after T ~ 1.8, the family /2 reproduces the exact 
result by itself, except for the vicinity of T ~ 2.7. As the family /3 gives relevant contributions 
in this interval, we plot the combined contributions of /2 and /3 in Fig. E^c). The resulting 
curve is clearly inaccurate and the reason is, once again, the proximity of these trajectories to a 
caustic, as can be seen from Fig. HJ^f). Thus, we use again the uniform formula Kf^{z"* , z' ,T) 
(now for families /2 and /3), and the resulting curves are shown in Fig. [2^c). Path C3 should 
be chosen for 1.8 < T < 2.55, and Ci, for T > 2.55, so that the exact curve is well reproduced. 

Fig. [2]^d) shows the uniform approximation for the whole time interval < T < 6, obtained 
by joining four pieces of the propagator: (a) families fl and /2 with path Ci for T < 0.65; (b) 
families fl and /2 with path C3 for 0.65 < T < 1.8; (c) families /2 and /3 with path C3 for 
1.8 < T < 2.55; (d) families /2 and /3 with path Ci for T > 2.55. 

We note that the trajectories of the /3 family do not seem to contribute to Eq. (jl6p . since 
their inclusion makes the results worse everywhere. However, as they become partners of the 
contributing trajectories belonging to /2 (when these approach the caustic), their contributions 
become fundamental to control the divergence of the family /2 when Eq. (|32p is used. 

7. The Nelson potential 

In the previous sections, both the theory and the numerical example were presented for 
one-dimensional systems only. Here we show an application for a two-dimensional chaotic 
system. The derivation of expressions equivalent to ()16p and (j32p for this case can be found in 
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Figure 3. Comparison between the exact prop- 
agator and the two semiclassical approaches for 
the Nelson potential. Panel (a) shows the exact 
curve (full line) and the individual contribution 
of two families of complex trajectories (dashed 

(2) 

lines for /I and dotted lines for /2) to Ksc . 
Panel (b) shows a blow up in the caustic region. 
The dash-dotted line corresponds to the com- 
bined contribution of the two families of trajec- 
tories. The uniform propagator is also shown 
for the paths Ci (squares), C2 (circles) and C3 
(triangles). Panel (c) shows the exact (full line) 
and final semiclassical results (dashed line). 



references |21l [24} [34] . The second order and uniform semiclassical formulas are given by 
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For the numerical application we have chosen the Nelson Hamiltonian, 

H = -{pl+pl) + {y-x^/2) +0.05x2. 

This system has been widely studied both classically [35j and quantum mechanically [21j. In 
particular, the coherent state propagator and its semiclassical quadratic approximation were 
investigated in great detail in [21] and some strongly divergent semiclassical behavior due to the 
presence of phase space caustics were identified. Here we revisit this problem, using the same 
parameters for which the caustics were found and apply the uniform formula. 

The widths of |z') and |z") were fixed to = = 0.2 and h = 0.05. The eight remaining 
parameters fixing the initial and final coherent states are: x' = x" = 0.72, y' = y" = 0.24, 
p'x = Px = —0.75 and p'y = Py = —0.63, which refer to a diagonal element of the propagator. 

Two families of trajectories (/I and /2) contribute to the propagator ([37|) when T is varied 
from 6 to 8.5, as shown in Fig. [3||a). Family fl reproduces very well the exact propagator 
for the range 6 < T < 7.2, while /2 does the same for 7.6 < T < 8.5. In the vicinity 
of T « 7.4, the semiclassical result shows a divergent behavior. As shown in Fig. ^h), the 
combined contributions of these families only makes things worse. A careful analysis of the 
classical phase space shows that this region is close to a caustic |21j . Therefore, it is indicated 
to use expression (|38p instead of (j37p . The results for the uniform approximation are shown in 
Fig. El^b) for the paths Ci, C2 and C3. Taking path Ci for T < 7.38 and C3 for T > 7.38 kills 
the divergence and the exact result is reproduced, as shown in Fig. [3|^c). 

8. Final Remarks 

Focal points and caustics are well known sources of inaccuracies in semiclassical approximations. 
The most systematic way to derive semiclassical formulas that avoid such divergences is to use 
the method of Maslov, which consists basically of two steps. The first step is to compute 
the semiclassical propagator in a dual representation. The caustics in the original and the 
dual representations lie generically in different regions of phase space, so that at least one of 
the propagators is well behaved at any given point. The second step is to transform back 
to the representation of interest doing the corresponding integral by the stationary phase 
approximation, but expanding beyond the second order. 

Coherent states lack a natural dual representation and we have defined one such represen- 
tation in [23]. Although the transformation leading from the Bargmann to the dual form is 
not a simple Fourier transform like in the case of position and momentum, it leads to a similar 
Legendre transformation of the action in the semiclassical limit. With this representation in 
hand the second step of the Maslov method can be carried out as usual. It is not clear at this 
point if this alternative representation can be useful in other contexts and work in this direction 
is in progress. 

Finally we remark that one of the few semiclassical approaches that naturally avoids caustics 
is the initial value representation of Herman and Kluk \36\ [371 EH] (see also [39]), for which 
the tangent matrix elements, that go to zero at the caustic, appear in the numerator and do 
not lead to divergencies. However, convergence problems due to highly oscillatory contributions 
have been reported for chaotic systems |40j . which also required the development of additional 
techniques and methods. 
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